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We study single- and two-atom van der Waals interactions of ground-state atoms which are both 
polarizable and paramagnetizable in the presence of magnetoelectric bodies within the framework 
of macroscopic quantum electrodynamics. Starting from an interaction Hamiltonian that includes 
particle spins, we use leading-order perturbation theory to express the van der Waals potentials 
in terms of the polarizability and magnetizability of the atom(s). To allow for atoms embedded 
in media, we also include local-field corrections via the real-cavity model. The general theory is 
applied to the potential of a single atom near a half space and that of two atoms embedded in a 
bulk medium or placed near a sphere, respectively. 

PACS numbers: 34.35.-|-a, 34.20.-b 42.50.Nn 



I. INTRODUCTION 

The dispersion interaction between neutral and (unpo- 
larized) atoms or molecules — commonly known as the 
van der Waals (vdW) interaction — is, together with 
Casimir-Polder and Casimir forces, one of the conse- 
quences of zero-point fluctuations in quantum electrody- 
namics (QED) (for a recent review, see Ref. [IJ). The in- 
teraction potential of two polarizable atoms in free space 
was first studied for small distances (nonretarded limit) 
by London using second-order perturbation theory 
In this limit the result is an attractive potential propor- 
tional to with r being the interatomic distance. The 
London formula was extended to arbitrary distances by 
Casimir and Polder using fourth-order perturbation the- 
ory within the framework of QED [3]. They found an 
attractive potential proportional to r^^ for large sep- 
arations (retarded limit) where the potential is due to 
the ground-state fluctuations of both the atomic dipolc 
moments and the electromagnetic far field. Casimir and 
Polder also considered the potential of a polarizable atom 
in the presence of a perfectly conducting wall '3]. The 
result is an attractive potential which shows a depen- 
dence in the nonretarded limit and is proportional to z^** 
in the retarded limit (z being the atom- wall separation) . 

In the three-atom case, a nonadditive term prevents 
the potential from just being the sum of three pairwise 
contributions; it was first calculated in the nonretarded 
limit [J2J5I and extended to arbitrary interatomic dis- 
tances [6|. Later, a general formula for the nonadditive 
TV-atom vdW potential was obtained by summing the 
responses of each atom to the electromagnetic field of 
produced by the other atoms 0] , or alternatively, by cal- 
culating the zero-point energy difference of the electro- 
magnetic field with and without the atoms Q. 

The theory was first extended to magnetic atoms by 
Feinberg and Sucher [9| who studied the retarded inter- 
action of two electromagnetic atoms based on a calcula- 



tion of photon scattering amplitudes; their results were 
later reproduced by Boycr [10] using a zero-point energy 
technique. In this limit, the interaction potential of a 
polarizable atom and a magnetizable one was found to 
be repulsive and proportional to r~^. Later on, Feinberg 
and Sucher extended their formula to arbitrary distances 
[ill] . In particular, in the nonretarded limit the potential 
of the mentioned atoms is found to be repulsive and pro- 
portional to r~'*. The retarded Feinberg-Sucher poten- 
tial was extended to atoms with crossed electric-magnetic 
polarizabilities on the basis of a duality argument p^ . 
For the single-atom case, the atom-wall potential, calcu- 
lated in Ref. 3J in the retarded limit, has been gener- 
alized to atoms with both electric and magnetic polar- 
izabilities [13], showing that a magnetically polarizable 
atom in a distance I from a conducting wall is repelled 
by that wall due to a potential proportional to A full 
QED treatment has been invoked to study the potential 
of an excited magnetic atom placed inside a planar cavity 
for all distance regimes [isj . 

In order to extend the theory of atom-atom interac- 
tions to the case of magnetoelectrics being present, the ef- 
fect of the bodies on the fluctuating electromagnetic field 
must be taken into account. A general formula expressing 
the vdW potential between two polarizable ground-state 
atoms in the presence of electric bodies in terms of the 
Green tensor of the body-assisted electromagnetic field 
was first obtained using linear response theory |14| and 
later reproduced by treating the effect of the bodies semi- 
classically [15j]- Recently, an analogous formula for two 
polarizable atoms interacting in the presence of magneto- 
electric bodies was derived by using fourth-order pertur- 
bation theory within the framework of macroscopic QED 
[T6| . it was later generalized to N atoms [l3|- For atoms 
that are embedded in a host body or medium, the local 
electromagnetic field experienced by them differs from 
the macroscopic one. Hence, the theory of vdW interac- 
tions must be modified by taking local-field corrections 
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into account. One approach to this problem is the real- 
cavity model, where one assumes that each guest atom 
is surrounded by a small, empty, spherical cavity p^ . 
It has been used to study local-field corrections to the 
spontaneous decay rate of an atom embedded in an ar- 
rangement of magnetoelectric bodies and/or media [l^ 
and was recently applied to obtain local-field corrected 
formulas for one-atom and two-atom vdW potentials of 
polarizable atoms within such geometries [20|. 

In this article, we generalize the theory of ground-state 
single- and two-atom vdW potentials in the presence of 
arbitrarily shaped magnetoelectric bodies to atoms ex- 
hibiting both polarizabilities and (para-) magnetizabil- 
ities. Such a theory includes and generalizes the re- 
cently studied potential of two polarizable and magnetiz- 
able bodies embedded in a bulk magnetoelectric medium 
pll |. This article is organized as follows. In Sec. Ull the 
multipolar atom-field interaction Hamiltonian is derived 
for atoms that are both electric and (para)magnetic. In 
Sec. Illli general expressions for single- and two-atom po- 
tentials are derived using perturbation theory. Local-field 
corrections are considered in Sec. IIVI while in Sec. |Vl we 
apply our theory by studying the examples of (i) an atom 
in the presence of a half space, (ii) two atoms in bulk me- 
dia, and (iii) two atoms in the presence of a sphere. A 
summary is given in Sec. IVII 



II. ATOM-FIELD INTERACTIONS IN THE 
PRESENCE OF SPINS 

The interaction of individual atoms with medium- 
assisted electromagnetic fields has extensively been dis- 
cussed for spinless atoms [l^, [13, H^l . In order to cor- 
rectly describe the paramagnetic properties of an atom, 
it is crucial to include the spins of its constituents in 
the considerations. A neutral atom (or molecule) A thus 
has to be regarded as being a collection of (nonrelativis- 
tic) particles a € A which in addition to their charges 
la iJ2a&Ala = 0)j uiasses iTia, positious , canouically 
conjugate momenta have spins s^. The particle spins 
give rise to magnetic dipole moments 7qSq, where is 
the gyromagnetic ratio of particle a ['je = —ege/(2me) 
for electrons with — e: electron charge; — 2, electron 
g-factor; mg: electron mass]. While leaving the atomic 
charge density 



PA 



(r) = ^ qaSir - r^) 



(1) 



and polarization 



generalize to 

jAir) = X! Y [f'aSir - ra) + 5{r - ra)ra 

a^A 

laSa X V5{r - To) 



(3) 



and 



a^A 



MA{r) = X! Y / ^'^'^ 5{r -TA- ctt-q) x r 



aeA 



(4) 



In Eqs. ([2]) and (|4]), Tq, — denotes the position of 
the ath particle relative to the center-of-mass position 



Erua . 
r 



aeA 



WLA 



(5) 



{mA = X^aeA "^a)' with the associated momenta being 



Pa = Pa 



TUA 



-PA 



(6) 



and Pa = Sas/iPa' respectively. Since the current den- 
sity associated with the spins is transverse, the continuity 
equation pa + ^ ■ Ja = remains valid. In addition, the 
atomic charge and current densities can still be related 
to the atomic polarization and magnetization via 



PA = - V • P4 , 
ja = Pa + V X Ma+ Jro 



(7) 
(8) 



as in the case of spinless particles, since the particle spins 
lead to equal contributions on the left and right hand 
sides of Eq. ([5]), as an inspection of Eqs. ^ and Q 
shows. In Eq. ([H]), 



3Ro 



i V X 



-P4 X f A - X Pa 



(9) 



is the Rontgen current density associated with the center- 
of-mass motion of the atom [2^ [2^. Further atomic 
quantities of interest are the atomic electric and mag- 
netic dipole moments 



^A = ^ QaTa = ^ qa-Ta 



aeA 



aeA 



and 



PA{r)=^qaral da 5{r - t-a - (rr^,) (2) rhA ^ y^\'^Ta x^a + laSa 

aeA ^0 aeA ^ ^ 



(10) 



(11) 



unaffected, the spin magnetic momentsdo contribute to 
the atomic current density and mag netization, so that the 
expressions given in Refs. [H, [2J| for spinless particles 



which emerge from the atomic polarization ^ and mag- 
netization ([3]) in the long-wavelength approximation, as 
we will see later on. The first and second terms in 
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Eq. (|lip obviously represent the orbital angular momen- 
tum and spin contributions to the magnetic dipole mo- 
ment. 

In order to account for the interaction of the spins 
with the magnetic field, a Pauli term has to be included 
in the minimal-coupling Hamiltonian given in Ref. [l6l | 
for spinless atoms interacting with the quantized elec- 
tromagnetic field in the presence of linearly responding 
magnetoelectric bodies, viz., 



duhhj fl{r,uj) ■ f\{r,uj) 



E 

aeuA 



E 



" Q./3euA ^' " 



(12) 



The first term in Eq. (|T2|) is the energy of the electromag- 
netic field and the bodies, expressed in terms of bosonic 
(collective) variables fx{r,uj) and f\{r,uj) (A,A'£{e,TO}, 
with e, TO denoting electric and magnetic excitations), the 
second term is the kinetic energy of the charged particles 
constituting the atoms, the third and fourth terms de- 
note their mutual and body-assisted Coulomb potentials, 
respectively, and the last term is the newly introduced 
Pauli interaction of the particle spins with the body- 
assisted magnetic field. Note that the scalar potential 
(p, the vector potential A, and the induction field B are 
thought of as being expressed in terms of the fundamen- 
tal bosonic fields fx and /j^ [lllil- 

To verify the consistency of the Hamiltonian (fT2|) . we 
need to show that it leads to the correct Maxwell equa- 
tions for the electromagnetic field and the Newton equa- 
tions for the particles. As in the case of spinless particles, 
the total electromagnetic field can be given by 



(13) 
(14) 



where E, B, D and H are the body-assisted electro- 
magnetic-field strengths [H, 113] , and 



w = E 



Qa 



aeA 



(15) 



is the Coulomb potential due to atom A. Since the atomic 
charge density ([T]) is not affected by the particle spins 
either, the Maxwell equations 



V • B = 0, 

A 



(16) 
(17) 



which are not governed by the system Hamiltonian, are 
not changed by the presence of spins. It is obvious that 
the Maxwell equation 



V X S + B ^ 



(18) 



also remains unchanged, because the Pauli interaction 
term commutes with the B-field and hence its inclusion 
does not lead to an additional contribution in Heisen- 



berg's equation of motion B— {\/h) 
by the commutation relation 27] 



H,B 



Mr), Mr') =miAr-r') 



As implied 



(19) 



[5^j{r), transverse delta function], the Pauli interaction 
does lead to an additional contribution 



^ -iaSa ■ B{ra),D{r) 

. aeuA 



aeuA 



in the Heisenberg equation of motion T> = {i/h) [H, X>] , 

which coincides with the spin-induced component of 
Jsecond term in Eq. ([3])]. Hence, the Maxwell equation 



(21) 



holds in the presence of spin when using the amended 
atomic current density 

Next, consider the equations of motion for the charged 
particles. Using the Hamiltonian IT^ . we have 



1 

TOn 



Pa - qaA{fa) 



(22) 



as in the absence of spins. Equation ((22)) implies that the 
Pauli interaction gives rise to a contribution 



- ^ -ypsp ■ B{ff3),mara 

/3GUA 



= 7a Vq 



Sa ■ Bifa) 



(23) 

{\/h) H , rriara ■ Combining this with 
the contributions from the spin-independent part of the 
Hamiltonian 12711 . we arrive at 



Sa ■ B{ra) 



(24) 



The first two terms on the right-hand side of this equation 
represent the Lorentz force on the charged particles while 
the third term is the Zeeman force resulting from the ac- 
tion of the magnetic field on the particle spins. We have 
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thus successfully established a Hamiltonian [Eq. (fT2|) ] de- 
scribing the interaction of one or more atoms with the 
electromagnetic field in the presence of magnetoelectric 
bodies which generates the correct Maxwell equations for 
the fields and the correct Newton equations for the par- 
ticles. 

Due to the rather large number of atom-field and even 
atom-atom interaction terms, the Hamiltonian (|12p may 
be not very practical as a starting point for calculations. 
As an alternative, we use the multipolar-coupling Hamil- 
tonian which for neutral atoms follows from a Power- 
Zienau-WooUey transformation [2^ 



6' = UOW 



with 



U = exp 



i/dV5:^.A 



(25) 



(26) 



upon expressing the Hamiltonian ()12|) in terms of the 
transformed variables. The only difference with respect 
to the case of spinless particles is the Pauli interaction 
term, which is invariant under the transformation since 
B' = B and s'^ = Sa- The multipolar-coupling Hamilto- 
nian can thus be given in the form of 



H 



(27) 



where 



E 

A— e,m 



dV / d^S^/^(r,c.)./;(r,c.), (28) 



4- 2to„ 2eo J 



2mA 



Pa 

2mA 



J2Er\n'A){n'A\, 



(29) 



with and E^' denoting, respectively, the eigenstates 
and eigenvalues of i?^, and 



H'af = - j d-V (p'a-E' + M'a-B') 

[ d?r^ -P'axB'+Y \ [ dVS' X B' 

J ruA 2m„ J 



aeA 



(30) 



with 



aeA 

-f^xf'J {r-r'A-<j¥'^)] + J2 iMr-f'^) (31) 



being the canonical magnetization and 



Kir) = P'Air) + g„r'„ / da a5 {v-v'a- ar'^) 

^A Jo 

J2 10^0 f da a5 {r-r'^- aT'p) . (32) 



m, 
mA 



In the long-wavelength approximation, the atom-field 
coupling Hamiltonian reduces to 



H'af = -d'A-E'{f'A)-rh'A-B'irU) + ^-d'AxB'if'A) 

ruA 



E: 

aeA 



K X B'ir'A) 



8mA 



d'A X B'ir'^ 



where 



and 



(I'a 



a£A 



aeA 



rriA = 



E 



Qa —I — / , ^/ 



(33) 



(34) 



(35) 



are the atomic electric and (canonical) magnetic dipole 
moments, respectively. In Eq. the first and the sec- 
ond terms are the electric and magnetic dipole interac- 
tion, respectively, the third term is the Rontgen interac- 
tion associated with center-of-mass motion, and the last 
two terms are generalized diamagnetic interactions. The 
Rontgen interaction becomes important when studying 
dissipative forces such as quantum friction [30l |. In this 
work, we are mainly interested in the interaction of atoms 
at given center-of-mass positions va featuring electric as 
well as paramagnetic properties which, upon discarding 
the last three terms in Eq. ([33]), can be described by the 
interaction Hamiltonian 



H'af = -d'A ■ E'irA) - rh'^ ■ B'{rA). 



(36) 



We conclude this section by recalling some relations 
that will be needed for the calculation of vdW potentials. 
The electromagnetic fields E' and B' are expressed in 
terms of the fundamental bosonic fields 



/;(r, uo)J2. ir\ ^')\ = SxySiu; - J)8{r - r'), 
/;(r,c.),A,(r',c.')] =0 
(A, A' G {e, to}) according to 

E'{r)= J2 y'dV^"d..GA(r,r',c^). A(r', 



(37) 
(38) 



OJ 



aeA 



H.C., (39) 



5 



/;(r',c^) + H.c., (40) 



where the quantities G\{r,r',uj) are related to the clas- 
sical Green tensor G(r, r', as 



Ge(r,r',w) = i^ 

c V '"■£0 



Im e(r',w)G(r,r',w), (41) 



(42) 

For an arbitrary arrangement of linearly responding mag- 
neto-electric bodies described by a permittivity e{r,Lu) 
and a permeability fi{r,Lu), the Green tensor obeys the 
differential equation 



„ 1 ^ , X 

Vx — -Vx - — e{r,uj) 



G(r, r', Lj) = J(i — r'), 
(43) 



has the useful properties 

G*{r,r',uj) = G{r,r',-uj*), (44) 
G(r,r',c^) = GT(r',r,o.), (45) 

and satisfies the integral relation [23| 



^ /'d3sGA(r,s,c.).Gf(r',s, 

\ 



TT 

The ground state |{0'}) of Hp is defined by the relation 
fx{r, w)|{0'}) = for all A, r, lj. Since we will exclusively 
work with the multipolar-coupling Hamiltonian, we will 
henceforth drop the primes indicating the Power-Zienau- 
WooUey transformation. 



III. VAN-DER-WAALS POTENTIALS 

According to the well-known concept of Casimir and 
Polder [31], vdW forces on ground-state atoms can be de- 
rived from the associated vdW potentials, which in turn 
can be deduced from a perturbative calculation of the 
position-dependent parts of energy shift induced by the 
atom-field coupling. 



A. Single-atom potential 

Let us consider a neutral ground-state atom A at a 
position va in the presence of arbitrarily shaped magne- 
toelectric bodies. With the atom-field interaction Hamil- 
tonian given by Eq. p6p (recall that we have dropped all 



primes) , the vdW potential of the atom follows from the 
second-order energy shift 



AE 



E 

/#0 



{0\Haf\I){I\Haf\0) 

El — Eq 



(47) 



where |0) = |Oyi)|{0}) denotes the quantum state where 
both atoms and the body-assisted electromagnetic field 
are in their ground states. Note that the summation 
in Eq. (|47p includes position and frequncy integrals. Re- 
calling the interaction Hamiltonian p6p , we see that only 
intermediate states |/) in which the atom is in an excited 
state and a single quantum of the fundamental fields is 
excited contribute to the sum and hence, Eq. (|47p may 
be specified as 



AE : 



k A= 



d-^r 



duj 



+ 



X \{QA\m\HAF\lx{r,uj))\kA) (48) 

[c.^ = iE\ - E\)/h, \lx{r,u)) = /|(r,c.)|{0})]. Using 
the expansions (I39p and ((40|) as well as the commuta- 
tion relations ([37|) and ([38|) . the matrix elements of the 
interaction Hamiltonian (1361) are found to be 



m^A ■ X Gx{rA,r,uj) 



{QA\m\HAF\lx{r.uj))\kA) 
^-df-G.irA^r^u) "'^ 



lOJ 



(49) 



[<= {OA\dA\kA), mf ^ {OA\niA\kA)]. 

With AE being quadratic in the matrix elements, there 
are three classes of contributions to the energy shift. The 
contribution involving two electric-dipole transitions is 
known to lead to the electric single-atom vdW potential 

M 



UeivA) 



2tt 



poo 

/ def tr[a^(iO-G«(r^,r^,iO] 
Jo 



= ^ J^"° d^eaAiiOtrG^'HrA^VAAO. (50) 
where G^^-* is the scattering part of the Green tensor, and 

9 , ,fe jOkjkO 



uj^ — iwe 



2 ^ clldO/p 



e^o 3^ ^ 



3h ^ (i^^a)^ — — iwe 



\ = aA{^)\ (51) 



(I denotin the unit tensor) is the atomic ground-state 
polarizability. The second lines in Eqs. ([50]) and ([5T|) are 
valid for isotropic atoms. 

The contribution AE„^ to AE which involves two 
magnetic-dipole transitions can be calculated by substi- 
tuting the second term in Eq. P^)) into Eq. (I48p and using 
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the integral relation (|46)) . resulting in 



duj 



■ V X Im G(r, r', w) x V' • m 



(52) 



[(T X V)^ = — VxT^]. The relevant, position-dependent 
part of Ai^m is obtained by replacing the Green ten- 
sor with its scattering part. After writing Im G = 
(G — G*)/(2i), making use of Eq. ((44|) . and transform- 
ing the integral along the real axis into ones along the 
purely imaginary axis (cf. Ref. [23j), the resulting mag- 
netic single-atom potential reads 



2tt 



d^ tr 



^^°°de/3^(iC)trL(i)(r^,r^,iC), 



where 



L(r,r',w) = V X G{r,r',bj) x V 



(53) 



(54) 



[note that L^^^ refers to G^^'], and 



2 X ^ 

/3a (w) = lim - 



■\^I3a{uj)\ (55) 



is the atomic ground-state magnetizability. The second 
lines in Eqs. ([55]) and ([55]) are again valid for isotropic 
atoms. 

We restrict our considerations to non-chiral atoms and 
molecules whose energy eigenstates can be chosen to be 
eigenstates of the parity operator. Contributions to the 
energy shift that contain one electric-dipole transition 
and one magnetic-dipole transition can then be excluded, 
since (Ia is odd and rriA is even under spatial reflection. 
Hence, the total vdW potential of a single ground-state 
atom that is both polarizable and (para)magnetizable 
and is placed within an arbitrary environment of mag- 
netoelectric bodies reads 



U{rA)^Ue{rA) + U,n{rA). 



(56) 



with Uf, and C/„i being given by Eqs. ([50)1 and re- 
spectively. To our knowledge, the magnetic part of this 
potential has been derived for the flrst time in this gen- 
eral form. 



B. Two-atom potential 

We now consider two neutral ground-state atoms A 
and B at given positions ta and rs in the presence of 



arbitrarily shaped magnetodielectric bodies. The two- 
atom vdW potential follows from the fourth-order energy 
shift 



{0\HAF + HBF\in){in\HAF + HBF\n) 



{Em — Eq) 



IJIJII^O 

^ {II\Haf+Hbf\I){I\Haf + Hbf\0) 
{Eji — Eo){Ej — Ef)) 



, (57) 



where |0) = jO^) \0b) |{0}) is the ground-state of the com- 
bined atom-fleld system. With the interaction Hamilto- 
nian being given by Eq. ([36]) . the summands in Eq. ((ST]) 
vanish unless the intermediate states |/) and \III) are 
such that one of the atoms and a single quantum of the 
fundamental fields are excited. The intermediate states 
\II) correspond to one of the following three types of 
states: (i) both atoms are in the ground state and two 
field quanta are excited, (ii) both atoms are excited and 
no field quantum is excited, and (iii) both atoms are ex- 
cited and two field quanta are excited. All the possi- 
ble intermediate states together with the respective en- 
ergy denominators are listed in Tab. |l] in App. |Al In 
addition to the matrix element an evaluation of 

Eq. (j57[) hence requires matrix elements of the interaction 
Hamiltonian (|36p which involve transitions of the body- 
assisted field between single- and two-quantum excited 
states. Recalling the definitions ((39|) and ([40]) as well as 
the commutation relations (13711 and (l38l). one finds 



{kA\{h^^^ (ri , wi) jiJ^F Ma^js {r2, UJ2) h^t, (ra, W3)) |0a) 



V2 



[df ■GxJrA,r2,L02)], 



'(12) r,ko 



V2 



[df-Gx,{rA,r3,u;,)]^ 



U!3 



[m'X"-VAxGxArA,r2,LU2)], 
[mf ■ Va X GA3(rA,r3,a;3)], 



W2V2 
ii5(i2) 



V2' 



(58) 



with \l,{r,Lu)ly{r',iu')) ^ ^flir,u;)fl,ir',co')\{0}) 
and 



^iiiu) = ^Xf^xji^t^ {rf, - rt,)6{ujf, - lu^). (59) 



The two-atom potential follows from those contribu- 
tions to the energy shift (|57p in which each atom under- 
goes exactly two transitions. As in the single-atom case, 
we distinguish different classes of contributions according 
to the electric or magnetic nature of those transitions. 
Those involving only electric transitions of both atoms 
are known to lead to the electric-electric vdW potential 
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X tr [a A (iC) • G (r^ , rs , ' (iC) ■ G (r^ , , i^)] 

X tr[G{rA,rB,iO ' G(ri3,rA,iO] (60) 

[recall Eq. ([5T|) ]. where the second equality is valid for 
isotropic atoms. 

Next, we calculate the electric-magnetic vdW poten- 
tial Uem, which is due to contributions of atom A under- 
going electric transitions and atom B undergoing mag- 
netic transitions. Each of the possible intermediate-state 
combinations listed in Tab. |T] contributes to J7emi where 
we begin with the intermediate states of case (1). Substi- 
tuting the respective matrix elements from Eqs. (|49p and 
([58|) into Eq. ([57|) and using the integral relation (I46p . we 
find 



Xt'^ io \Dla Du 



(61) 



x| df ■lmG{rA,rB,oj) xVb -m^B 
X [m% ■ Vb X ImG(rB,r^,c^') • df] 



with the energy denominators Dia and Dit being given in 
Tab. m Without loss of generality, we have assumed that 
the matrix elements of the electric- and magnetic-dipole 
operators are real-valued quantities. One can then easily 
find that the contributions AE^^ {k G {2,3,..., 10}) 
from the other possible intermediate-state combinations 
differ from Eq. (|61|) only with respect to their energy 
denominators and signs. Case (6) leads to two terms 
with different energy denominators l/Dga + iZ-Deb, just 
the same as case (1), while all other cases only give rise to 
a single term each. Furthermore, the contributions from 
cases (3)-(5), (8)-(10) differ in sign from Eq. The 
electric-magnetic vdW potential can be found as the sum 

of all contributions Uem [fA TfB) = J2k ^^em ■ In analogy 
to Ref. lial it can be seen that the denominator sum 



(62) 



1 1 


1 


1 


1 


1 






~Ds' 






1 1 


1 


1 

















can be replaced by 



(63) 

under the double frequency integral in Eq. (|6ip , where we 
have used the definitions of the denominators in Tab. U 



and exploited the fact that the remaning integrand is 
symmetric with respect to an exchange of to and uj' . This 
results in 



C/e™(r^,rs) = -|4E^ 



{^1 



Jo 

1 



do; 
1 



dJ 



[x))(a;^+cj) ^ u' UJ — uj' 



X {[df ■linG{rA,rB,oj) x Vij • m^'] 

X [m^. Vb xImG(rB,r^,L^')-<]}- (64) 

The integral over uj' can be performed by using the iden- 
tity ImG = (G - G*)/{2i) and Eq. to yield [l| 



1 



UJ + Oj' 



ImG(rB,rA,cj') 



= --uj[G{rB,rA,uj) + G*{rB,rA,uj)] ■ (65) 

After substituting Eq. into Eq. and transform- 
ing the w-integrals by means of contour-integral tech- 
niques to run along the positive imaginary axis, one ob- 
tains 

UUva^vb)^^ TdCf 
X tr \aA{iO ■ KT(rs,r^,ie) •/3s(iO • K(rs,r-^,ie) 

-^l^^d^eaAmPsiio 

X tr K"^(rB,rA,i^) • K(rB,rA,i^) , (66) 



where 



K{r,r',uj) = V X G(r,r',cj), 



(67) 



and the second equality holds for isotropic atoms. Ob- 
viously, the magnetic-electric potential J7me('"yi, f"B), 
which is due to all contributions of atom A undergo- 
ing magnetic transitions and atom B undergoing electric 
transitions, can be obtained from Eq. (|66p by interchang- 
ing A and B on the right hand side of this equation. The 
magnetic-magnetic potential Umm, associated with mag- 
netic transitions of both atoms, can be found in a pro- 
cedure analogous to the one outlined above for deriving 
Eq. ([5^ . resulting in 



Unini{rA,rB) 



h44 

X tr [I3a{M) ■ L(rA, rs, i^) • Pb^ ■ L(rB, r^, i^)] 

,2 



X tr [L(rA, Tij, if ) • L(rB, r^, i^] , (68) 
where the second equality again holds for isotropic atoms. 
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We have thus calculated all those contributions to the 
energy shift where both atoms undergo exactly two tran- 
sitions of the same type (electric/magnetic). The re- 
maining contributions of one or both atoms undergoing 
an electric and a magnetic transition can again be ex- 
cluded from a parity argument for the non-chiral atoms 
under consideration in this work (for the interaction of 
two chiral molecules in free space, see Ref. [33]). The 
total two-atom vdW potential of two polarizable and 
(para)magnetizable atoms placed within an arbitrary en- 
vironment of magnetoelectric bodies is hence given by 

U{rA,rB) = Uee{rA,rB) + Uern{rA,rB) 

+ U,ne{rA,rB) + Ujnm{rA,rB), (69) 

together with Eqs. ([60]) . ([66|) and (l68|) (the diamagnetic 
contribution to the dispersion potential of two atoms in 
free space is discussed in Refs. [ssl Is^. |35|V 



IV. LOCAL-FIELD CORRECTIONS 

The single- and two-atom potentials given in Sec. IIIII 
refer to atoms that are not embedded in media, i.e., 
£{rA(B),^) = A'('"A(B), ^) = 1- When considering guest 
atoms inside a host medium, one needs to include local- 
field corrections to account for the difference between the 
macroscopic electromagnetic field and the local field ex- 
perienced by the guest atoms. A possible way to treat 
local-field effects is offered by the real-cavity model jTsj l. 
where small spherical free-space cavities of radius Rc 
surrounding the atoms are introduced. As shown in 
Ref. [l^ , the local- field corrected forms of the Green ten- 
sor read, in leading order of uRc/c, 



Gioc{rA,rB,uj) 



Sea 



2eA + l 



G{rA,rB,uj) 



3£f 



2eB + l' 



(70) 



[we have discarded the position-independent first term 
on the right-hand side of Eq. ([7T|) ] and 



Uee{rA,rB) 



2vr Jo 
3eA(iO 



26^(10 + 1 



3£B(iO 



2es(iO + l 



X tr [a^ (iO • G (r A , rs , i^) • (i^) • G (r^ , r ^ , i^] 



(73) 



For magnetic atoms the vdW potentials depend on 
spatial derivatives of the Green tensor. Hence, the re- 
spective local-field corrected tensors cannot be derived 
directly from Eqs. ([70]) and ([7T|) . because the correction 
procedure does not commute with these derivatives. As 
shown in App. [B] the required local-field corrected forms 
of the tensors L [Eq. ([54])] and K [Eq. ^7}] within leading 
order of uiRc/c are given by 



3 ^^(5£A-1)-3^A-1 c 



UjRr, 



5 {2^iA + l)^ 

2 

L'^^\rA,rA,uj 



3/x^n 



1 



(2/iA + l)2 3 



^ (2MA+1 

Uoc{rA,rB,uj) 
and 



\-[rA,rB,uj) 



2flA + 1 



2/iB + 1 



I 

(74) 
(75) 



G,^,{rA,rA,u:)^ — 



2eA + l uj^R^ 



3 £^(5/1^-1 

'5 



-3e^— 1 c 



(2eA + l)^ 



ivRc 



3e^n' 



{2eA + l)^ 



+ 



3£a 

2SA + l 



G''-^\rA,rA,uj), 



(71) 



where eA{B) = e(^A(s),^) and fiA{B) = /^(^yi(s), w), re- 
spectively, are the permittivity and permeability of the 
unperturbed host medium at the position of the guest 
atom A{B) (nA{B} = ^/£AXB)pA(Bj) and G is the uncor- 
rected Green tensor. Inserting the corrected Green tensor 
into Eqs. ([50| and ([60|) . one obtains the local-field cor- 
rected electric contributions to the single- and two-atom 
vdW potentials ^ 



TT / \ ^/^O 
Ue{rA) = -7^ 
ZTT 



3£A(iO 



n 2 



tr[aA(ie)-G('^(r^,r^,iC)l (72) 



Kioc(ryi,fB,cj) = - — - K(rA,rB,w) ^^■^ 



2MA + 1 



Replacing in Eq. ([531) L'^^ with Lj^^^ from Eq. ([74]), we 
obtain the local-field corrected magnetic single-atom po- 
tential 



2£i 



(1) 



(76) 



2^A(iO + l. 
trf/3^(ie)-L(i)(r^,r^,id , (77) 



where a position-independent term has been discarded, 
as in the electric case. To obtain the local-field corrected 
contributions Uem and Umm to the two-atom vdW poten- 
tial, we replace K and L with Kioc and Lioc in Eqs. 
and respectively, leading to 
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2 /-co 



UUrA,rB)^^l dee 



2n 



3eA(ie) 



tr 



and 



(78) 



n 2 r 



2/XA(ie) + 1 



2^B(ie) + 1 



tr[/3^(iO • L(r^,rs,iC) -/SsaO ' L(rB,r^,ie)] . (79) 



Recall that C/me('"A,^B) can be obtained from Eq. ([75)1 
by interchanging A and B on the right-hand side of this 
equation. Needless to say that Eqs. dTT]), (|75)) . 

and dZHl) reduce to Eqs. ([501), §Ui), ([53]), (El]), and 
respectively, when the atoms are situated in free space so 
that £a(B) = fJ-A{B) = 1- 



V. EXAMPLES 

We now apply the theory to some illustrative exam- 
ples and compare the results with the familiar results for 
nonmagnetic atoms, with special emphasis on whether 
the total potentials for electromagnetic atoms are invari- 
ant under a global duality transformation c^a^ (3 
[sl] . It will turn out that atoms situated in free space do 
respect this symmetry for the examples studied, while 
atoms embedded in media only do when the local-field 
corrections are taken into account. 



A. Single-atom potential: Half space 

First, we consider an isotropic atom A at a distance za 
away from a magnetoelectric half space of permittivity 
e(w) and permeability /J.(w) and choose the coordinate 
system such that the z-axis is perpendicular to the half 
space that occupies the region z < 0. Assuming that r 
and r' refer to two points in the free-space region z > 0, 
we have [33 



"/Lt(w)6- bo 



i(w, -r — w 

e ^ + 



8n% 



e{Lu)b-bo 



jiiu)b + bo e{Lo)b + bo P P 

[w^ =q±ibez, g-Le^], where 



P-^, bo^ ^q^-n^iiv)^ 



(80) 



(81) 



[q — \q\, n{uj) — ^/e(}^j)JI{Lo) , Re&, Re&o > 0], and the 
polarization vectors and are defined by (sq^q/q) 



OJ 



-(qe^ T i^eg 



(82) 



As shown in Ref. [sj, substitution of G^^-* from 
Eq. (|80|) into Eq. (|50|) yields for the electric part Ue of 
the single-atom vdW potential 



poo nOO 

Ue{rA) = 'i^ d^e^AiiO dq 
Jo Jo 

'fi{i^)b-bo e{i£_)b-bo 



1 f,-2bzA 

b 



fi{iOb + bo e{iOb + bo 



1 



2? |2 



(83) 



In the nonretarded limit of the atom-surface separation 
being small with respect to the characteristic atomic and 
medium wavelengths, Eq. (|83p simplifies to 



UeizA) - 

^J'ofl 



deaA(ie) 



IGtt^za 



deraA(ie) 



+ 1 



£(iO + 1 /i(ie) + 1 

2£(iO[n2(ie)-l] 



m) + 1? 



(84) 



In contrast, in the retarded limit of large atom-surface 
separation one finds that 



UeizA) 



3hcaA{0) 



dv 



1 



647r2eoz^ 

e(0)w- V"^(0)-1-Hw2 I ;x(0)w-^n2(0)-l-|-7 



£(0)w-f A/n2(0)-l + i;2 wV(0)t^+\A?(OPT+^ 

(85) 

To calculate the magnetic part C/„i of the single-atom 
vdW potential, we first combine Eqs. ([50]) and ([51)1 to 



d^q- 



e{uj)b - bo 
e{uj)b + bo 



e.es 



S-K^b 
H{uj)b- bo + _ 
^i{uj)b + bo ^P 



(86) 



Comparing Eqs. ^ [together with Eq. jMI)] and ^ 
[together with Eq. ((80|) ]. we see that the magnetic part 
Um can be found from the electric part Ue in Eq. ([55]) 
by replacing a a and e, with (3a/c^ and /i, respectively, 
in agreement with the duality principle [36[ . Needless to 
say that this symmetry also holds for the retarded and 
nonretarded limits. 
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B. Two-atom potential: Bulk medium 



and 



As a second example, we consider two isotropic atoms 
A and B embedded in an infinitely extended bulk 
medium of permittivity and permeability ^(w). To 
illustrate the relevance of the local- field corrections, let 
us first consider the uncorrected two-atom potential. By 
using the bulk-material tensors as given in Eqs. (|B3|1 and 
(|B7p . and calculating 



(1 — ikl)ei X I 



(87) 



{I = tb — ta, I = 6/ = Z/Z), which follows from 
Eq. dnZD together with Eq. dBSj, the potentials pi]) . 
and (1681) take the form 



UeeirA,rB) 

h 



167r3e2;6 



9[nim/c] 



Uem{rA,rB) = 

and 

Umm{rA,rB) 

where 



de^2^A(iO/3B(iO 

X (89) 



de/3^(iO/3s(iO 
X M'(ie)5["(iOeVc] , (90) 



c,(a;) =e"2"=(3-f 6x-|-5a;2-f 22;3-f x^), (91) 
/i(2;) =e-2'^(l + 2a; + a;2). (92) 

We see that due to the factors e^^(i^) and /i^(i^), the 
uncorrected quantities Uee and J7mm do not transform 
into one another under the duality transformation e <-!■ 
/i, c^a ^ (3. The same is true for the pair Uem and 
Ume- As a consequence, the uncorrected total two-atom 
potential ([55]) violates duality symmetry. 

By contrast, the local- field corrected two-atom poten- 
tial does obey the duality symmetry. From Eqs. ((73)l . 
([75)1 and dZH) [together with Eqs. jHS]), dHZ]) and (|B7)l ] 
we find that 



Uee{rA,rB) = - 



[2e(iO + 1] 



aA(iC)aB(iC) 

g[n(iOeVc], (93) 



UemirA,rB) = 



81e2(i^)^2(i^) 



da2^A(iO/?B(iO 



[2e(iO + lFMiO + l]2 



;i[n(iO^Vc], (94) 



Umni(rA,rB) 



167r3/6 



de/3A(ie)/3i3(iO 



^ pSiyi^^t"^^^^^^^^ (95) 



Inspection of Eqs. ((93l) - ((95| then reveals that the duality 
transformation £ ^ /i, c^a ^ /3 results in 



UeeirA,rB) ^ Umm{rA,rB), 
Uem{rA,rB) ^ Ujne{rA,rB), 



(96) 
(97) 



so the total vdW potential (l69|) is invariant under the 
duality transformation. The result clearly shows that (i) 
the inclusion of local-field effects is essential for obtain- 
ing duality-consistent results and that (ii) the real-cavity 
model is an appropriate tool for achieving this goal. 

It is instructive to inspect the nonretarded and re- 
tarded limits of Eqs. ((93)) - ([95]) . In the nonretarded 
limit where the atom-atom separation is small in com- 
parison to the characteristic atomic and medium wave- 
lengths, the approximations g[n{i^)^l/c] ~ g(0) and 
h[n{i£,)£,l/c\ ~ h{0) result in 



Uee{rA,rB) 



-3h 



dCaA(iOas(iC) 
8l£2(i^) 



Uem{rA,rB) 



(98) 



167r3Z4 7p 



[2e(iC) + l]2[2M(iO + l]2 



(99) 



167r3Z6 



U,nmirA,rB) = ^3 / d? /?A (iO/^B (if ) 

8lAi^(iO 



[2M(iO + 1]^ 



(100) 



In the retarded limit, the quantities a, /3, e, and /i can 
be replaced by their static values, leading to 



Uee{rA,rB) 



Uem{rA,rB) 



23ftcaA(0)as(0) 8l£2(0) 



64^3£2;7 ,^(0)[2e(0) + l]4 ' 

(101) 



7hcfioaA{0)PBiO) 



81n(0) 



Urnm{rA,rB) 



[2£(0) + l]2[2Ai(0) + l]2 

23hcfilPAiO)pB{0) 81^^2(0) 



, (102) 



64773^7 n(0)[2^(0) + 1]4 ■ 

(103) 
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Compared with two atoms in free space, one notices 
that the medium modifies the magnitudes of the inter- 
atomic potentials but does not change their signs. In- 
spection of Eqs. (|93p and ([M]) reveals that the medium 
always leads to a reduction of Uee and Umm- In the 
nonretarded limit, Uee is only influenced by the electric 
properties of the medium and Umm only by the magnetic 
ones [cf. Eqs. and (|100p ]. In contrast, Uem and 
Ume are diminished by the medium in the retarded limit, 
Eq. (|102p . but are enhanced by a factor of up to 81/16 
in the nonretarded limit [cf. Eq. 

In the retarded limit, the influence of the medium on 
all four types of potentials is very similar. The cou- 
pling of each atom to the field is screened by a factor 
9£(0)/[2£(0) -1-1]^ for polarizable atoms, and a factor 
9yLt(0)/[2/x(0) -|- l]'^ for magnetizable atoms. In addition, 
the reduced speed of light in the medium leads to a fur- 
ther reduction of the potential by a factor n(0). 

It should be pointed out that the uncorrected poten- 
tials Uem and Umm as given by Eqs. ([55]) and (^01) dif- 
fer from the corresponding results given in Ref. 21] by 
factors of and /i^^, respectively. The discrepancy 
is due to the different atom-field couplings employed: 
While our calculation is based on a magnetic coupling 
of the form m - B, a. m - H coupling is used in Ref. [2l[. 
The potentials derived therein thus do not follow from 
a Hamiltonian that is demonstrably consistent with the 
Maxwell equations and generates the correct equations 
of motion for the charged particles inside the atoms, 
whereas both of these requirements have been verified 
for the Hamiltonian (gT]) together with (gS]), and 
P6p employed in this work. Furthermore, in spite of the 
use oi&m- H coupling, the contribution due to the noise 
magnetization contained in H (cf. Ref. [l^H^]) was not 
discussed. The discrepancy would not have been notice- 
able if local-field corrections had been taken into account 
in Ref. [2l[ : When applying local-field corrections to the 
potentials stated therein, one recovers our local-field cor- 
rected Eqs. ([M|) and ([55)1 since the appropriate magnetic 
local-field correction factors are 3/x(i^)/[2/x(i^)-|-l] in that 
case, as opposed to the factors 3/[2/i(i^) -I- 1] arising in 
our calculation. 



C. Two-atom potential: Sphere 







z 




\ T 

\e 








\ y 


X 






FIG. 1: Two atoms A and B in the presence of a sphere 



U{rA, re) = f/^"^ (ta, vb) + U'-^\rA, tb) + U^^\rA, tb). 

(104) 

The potential contributions arising from the free-space 
part of the Green tensor can be found from Eqs. (|93p . 
(IM)) . and ([95)1 by setting e = ^ = 1. In the body- induced 
part of the interaction potential 

U\rA, tb) = U^^\rA, tb) + C/^^^ (r^, tb), (105) 

which arises from the scattering part of the Green tensor, 
the contributions uiV and t/ie'' to U^ can be taken from 



Ref. 1381 , and the contributions U, 



(1) 



and Umm to U^ can 



then be obtained from ui'^ and ui^^ by the transforma- 
tion a — !■ e ^ /i, as sketched in App. [Cl We may 
therefore focus on the calculation of the body-induced 
mixed contributions 



Ui^ivA^VB)^^ rdU'aAd 
Jo 



iO/3B(iO 



xtr K'-"^^{rB,rA,iO-^^"irB,rA,iO , (106) 



(1)/ 



Finally, let us consider two isotropic atoms A and B in 
the presence of a homogeneous sphere of radius R, per- 
mittivity e(Lu) and permeability /i(ci^). According to the 
decomposition of the Green tensor into a free-space part 
and a scattering part, each contribution to the two- atom 
vdW potential U{rA,rB), Eq. can be decomposed 

into three parts labeled by the superscripts (0), (1), and 
(2) , respectively, denoting the contribution from the free- 
space part of the Green tensor, the cross term of the free- 
space part and the scattering part of the Green tensor, 
and the scattering part of the Green tensor. 



Ui'2irA,rB)^^ raa'aAmM^O 

KWT(rB,r^,iO • K«(rB,r^,id • (107) 



X tr 



For this purpose, we choose the coordinate system such 
that its origin coincides with the center of the sphere 
(Fig. [T|). The scattering part of the tensor K{rB,rA,i^j) 
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can be given in the form (App. \C\\ 
ifco 



X sin 8 e.rj3 e^^ 



1 



^(2n + l) r^i?fQ„i^;(7) 



n{n + T) i^AB^l'Q^Fnh) ~ rsB^^ Qi 
X ^'n (7)] 69^60^4 +rBB^Q„P^(7)sinee0^e^^ 

n(n+l) '■ J f- J 

(108) 

[fco = w/c; rA(B) = |''A(B)|; 7 = cosO; e = 6*^ + 6*5, 
angular separation between the two atoms with respect 
to the origin of the coordinate system] , where 



M(^)[yojn(yo)]'jn(y) - [2/jn(?/)]'jn(yo) 



/i(w)[?/o/in^(yo)]'jn(y) - [yjn(2/)]'/in ^(yo 



(1)/ 



(109) 



£{^)[yojn{ya)]' jn{y) - [yjn{y)]' jn{yo) 



e{u;)[yohi'\yo)]'jn{y) - [yjn{y)]'h\," (yo) 



Qn^hlp{korA)h'^HkorB), 

Qf,^hip{korA)[zhllHz)l=kor,, 

Fn{x) = n{n + l)Pn{x) - xPi{x) 



(110) 

(111) 

(112) 
(113) 
(114) 

I 



[Pn[x), Legendre polynomial; yo = ^o-R; y = n{uj)ya\- 
Further, e^, eg, and are the mutually orthogonal unit 
vectors pointing in the directions of radial distance r, po- 
lar angle 6*, and azimuthal angle 4>, respectively (Fig. [1]). 
In order to facilitate further evaluations, it is convenient 
to represent the free-space part K'-^-', which can be ob- 
tained from Eq. (|87p for ^ = 1 and fc = feg, in the same 
spherical coordinate system as the scattering part. 



+ iBeeB^ipA + sin9 60^6^^ + lAe^^eg^), (115) 



where ^^(^b) is the component of I in the direction of 
^"a(-tb). 



Ia = fB cosQ — rAj Ib = tacosQ — rB- (116) 



Using Eqs. (fTU5)) and (fTTS]) in Eqs. (fTTO]) and pUT]) . we 
derive 



(2n+ 1) 



^^°"da3^^(iO/3B(iOe-'«/= (l + l) {n(n+l)sin2e[riBf (iC) 



n—l ■ 

r| 5^(10] Q„P;(7) + r^Sf (iOQf [^s^n(7) " ^^^^(7)] + rBB^i^Qn UFnin) - IbBM] , (117) 



Ui^l{rA,rB) 



^^M^ / dee^a^(ie)^B(iC)K(n' + l)n(n + l)Q„,g„sin^e 



' A' B „.-;^i + i)'^('i + 1) ^ 

X P:Al)Pn(.l) KB^' mBn'm + r|B,^(iOi3,r (iO] + (iOS^(iOOn'Q;t + ^i^^' (i^Sf (i^Qf'Qf ] 

X [^^„.(7)F„(7) + P:A7)Pnil)] - ^TaVbB^^ mBnimn'Qn K'h)Fnil) + (7)^n' (7)] } ■ (US) 

I 



As before, Ume{rA,rB) and Ume{rA, vb) can be obtained 
from Eqs. (|117p and (|118p by interchanging A and B. In- 
spection of Eqs. (I117|) and ()118p reveals that this is equiv- 
alent to the interchanging a P/A and e /i, which 
shows that the combination UemifA, ^b) + Umei'i'A, i'b) 
is invariant under the duality transformation. Recalling 
that UeeirA,T^B) + Umm{fA,fB) also obcys the duahty 



symmetry, the total potential U {ta, '"b) is duality invari- 
ant. 



Further analytical evaluation of the body-induced part 
of the potential is possible in the limiting cases of large 
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and small spheres. In the case of a large sphere, 

5a' = rA'-R^R, {A' = A, B) (119) 

/ < i? ^ e < 1 (120) 



[where the second condition in Eq. ()120|1 follows from 
the first one by virtue of 2i?sin(0/2) < I, cf. Fig. [T], we 
derive (App. ID|) 



Uinil'A.rB) 



{2/+(/+ + <5+)2 



X [{X^ - 6^5+) Jio + {X^ + 5-6+) Joi] + /3(2/^ + X^) 
(J20 + J02) + - /+J+) Jii} , (121) 



where X = RQ, 5±=5b±5a, 1+ = \ X'^ + Si, and 



Jkl = 



e(iO + 1 



LM(iC) + l. 

(122) 

In the case of a small sphere, R ^ rA'{A' = A,B), the 
main contribution to the frequency integrals in Eqs. (|117|) 
and (jllSp comes from the region where f ^ c/R, so that 
C/g^ can be approximated by the term n = 1 in Eq. (|117[) 
(cf. Ref. [H), leading to 



UemirA,rB) = 



A' B "'0 



da2„A(iO/3B(ie)e-('--*+'--+')«/^ (1 + IC/c) 



X I [2rB (1 + r^e/c) sin^ e + (^B cose) a(r^e/c)] {1 + rB^/c)rBa,p{iO 

+ [2rAil + TB^/c) sin^ O + iU- h cos 6) a{rB^/c)] (1 + rA^/c)rA^^ 

I 



, (123) 



where 



e{uj) + 2 



/3sp(w) 



47r 



/i(w) 



1 



i?'^ — -, — ^ ;r , 

Ho n(uj)+2 



and 



a{x) 



1 



(124) 
(125) 

(126) 



It is worth mentioning that the non-additive interac- 
tion potential of three atoms [polarizable atom A, mag- 
netizable atom B, and a third atom C of polarizability 
ac{uj) and magnetizability /^(^(w)] in free space may be 
obtained from Eq. (|123p by replacing asp(w)^ac(w) and 
Psp{uj)~^ Pc{c^)- By adding U'^^{rA,rB) from Ref. (H 
and U!^-i^j^{rA,rB) (cf. App. [Cjl . one can obtain the non- 
additive potential of three atoms, each being simultane- 
ously polarizable and magnetizable. 

Let us finally present some numerical results illustrat- 
ing the effect of a medium-sized magnetoelectric sphere 
on the vdW potential of two two-level atoms with equal 
transition frequencies. We again focus on the case where 
atom A is polarizable and atom B is magnetizable. The 
corresponding results for two polarizable atoms are given 
in Ref. [16] , from which, by duality, the analogous results 
for two magnetizable atoms can be inferred (see App. lCj) . 

Figures [Hand [3] show the ratio Uem/Uem obtained by nu- 
merical computation of Eq. (|104p together with Eqs. 



(for £ = ^ = 1), (HSl), (fTTT)) . and (fTTS)) . with the permit- 
tivity and permeability of the sphere being approximated 
by single-resonance Drude-Lorentz models, 



e{Lu) = 1 + 



H{uj) = 1 + 



(127) 
(128) 



In Fig. [21 two atoms at equal distances va = '"s from 
an electric sphere are considered and the ratio Uem/Uem 
is shown as a function of the angular separation & of 
the atoms, for three different values of the atom-sphere 
separation. It is seen that the presence of the sphere 
can lead to enhancement or reduction of the potential, 
depending upon B. To be more specific, Uem/Uem first 
increases with 0, attains a maximum, and then decreases 
with increasing S to eventually become minimal at = 
TT when the atoms are positioned at opposite sides of 
the sphere. Whereas the position of the maximum shifts 
with the atom-sphere separation, the minimum is always 
observed at = tt. Note that a magnetic instead of an 
electric sphere would lead to the same behaviour, because 
of duality. 

Figure [31 shows the dependence of the ratio Uem/Uem 
on the separation distance / between the two atoms for 
a configuration where the atoms are on a straight line 
through the center of a sphere (i.e., = 0), with the 
polarizable atom A being closer to the sphere than the 
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FIG. 2: The vdW potential of a polarizable and a magneti- 
zable two- level atom (transition frequency wio) in the pres- 
ence of an electric sphere of radius R = cjiova {ujp^/ujio = 3, 
ujt^/i-ou) ~ 1, 7e/i^io = 0.001) is shown as a function of the 
angular atom-atom separation Q. The values of ta ~ rs are 
1.03c/cJio (solid line), l.Sc/oJio (dashed line), and 2c/loiq 
(dotted line). 



magnetizable atom B. Note that in contrast to the previ- 
ous configuration, in this case an electric and a magnetic 
sphere do not lead to equivalent results by means of du- 
ality, because the positions of the electric and magnetic 
atoms are not equivalent. From Fig. [3lja) it is seen that 
in the case of an electric sphere the interaction potential 
is reduced compared to its value in free space; the ratio 
Uem/Uem dccreases vifith increasing I and approaches an 
asymptotic limit that depends to the distance between 
atom A and the sphere. In contrast, from Fig. EJ^b) it is 
seen that in the case of a magnetic sphere the interaction 
potential is enhanced compared to its value in free space, 
and a pronounced maximum of the ratio Uem/Uem is 
observed. For large atom-atom distances, Uem/Uem ap- 
proaches an asymptotic limit that is independent of the 
distance between atom A and the sphere. 



VI. SUMMARY AND CONCLUDING 
REMARKS 

We have extended the framework of macroscopic QED 
to paramagnetic atoms by introducing a Pauli term in 
the atom-field interaction. We have verified the consis- 
tency of our generalized Hamiltonian by showing that it 
generates Maxwell's equations and the correct equations 
of motion for charged particles with spin. On the basis 
of this Hamiltonian, we have employed leading-order per- 
turbation theory to generalize the theory of body-assisted 
one- and two-atom van der Waals potentials of polariz- 
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FIG. 3; The vdW interaction potential of two atoms with 
parameters as in Fig. [2]in the presence of (a) the same electric 
sphere as in Fig. [21 and (b) an analogous magnetic sphere is 
shown as a function of the atom-atom distance I for = 
and vb = ta + I- The values of ta are 1.03c/li;io (solid line), 
l.lc/tJio (dashed line), 1.3c/a;io (dotted line). 



able atoms to those that are both polarizable and mag- 
netizable. It is seen that, with respect to each atom, the 
generalized potential can be considered as a superposi- 
tion of contributions associated with the atomic polar- 
izabilities and magnetizabilities. We have extended the 
scope of our theory to atoms that are embedded in me- 
dia by implementing local-field corrections via the real- 
cavity model. We have found that local-field effects give 
rise to correction factors that depend on the permeability 
of the host medium for magnetizable atoms rather than 
the permittivity, as is the case for polarizable atoms. 

We have applied the theory to the single-atom poten- 
tial of an atom in the presence of a magnetoelectric half 
space and to the two-atom potential of atoms embed- 
ded in a bulk magnetoelectric medium or placed near 
a magnetoelectric sphere. The potential of a magnetiz- 
able atom in the presence of a half space has been found 
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to be very similar to the known respective potential of 
a polarizable one. We have shown that a bulk medium 
does not change the sign of the two-atom interaction, but 
can lead to enhancements and reductions, whereby in the 
nonretarded limit the potentials of two polarizable or two 
magnetizable atoms is only influenced by the electric and 
magnetic medium properties, respectively. For the two- 
atom potential in the presence of a sphere, the case of 
two magnetizable atoms was demonstrated to be analo- 
gous to the known case of two polarizable, so we have 
focussed on the sphere-assisted interaction of a polariz- 
able atom with a magnetizable one. We have obtained 
analytic results for a very large sphere (in which case 
the potential coincides with that of a half space) and a 
very small sphere (where the potential is analogous to 
the nonadditive three-atom interaction potential in free 
space, with the sphere taking the role of a third atom). 
Numerical results have been obtained for medium-sized 
spheres, where the sphere gives rise to enhancements and 
reductions of the potential, depending on the gemoetric 
arrangement of atoms and sphere: In particular, when 
the atoms are placed at equal distances from the sphere, 
the potential is enhanced (reduced) for small (large) sep- 
aration angles between the atoms, while a linear arrange- 
ment of the atoms and the sphere (with the polarizable 
atom being closer to the sphere) leads to reduction (en- 
hancement) for a electric (magnetic) sphere. 

For the examples involving atoms in free space, we 
have explicitly verified invariance with respect to a global 
interchange of e ^ /i and c^a /3, in a gree ment with 
the duality properties investigated in Ref. [3g| . The case 
of two atoms in a bulk medium has further revealed that 
this duality invariance only holds when accounting for 
local-field corrections. 



APPENDIX B: LOCAL-FIELD CORRECTED 
TENSORS L AND K 

The local-field corrected version of the tensor L de- 
fined by Eq. (|54p can be derived in complete analogy to 
the derivation of Eqs. ([70)1 and (|7ip . which was given in 
Refs. jl^, [20]. For this purpose we recall that the first 
term in Eq. (j7ip . i.e., the position-independent part of 
G[ll{rA,rA,Ltj), stems from the scattering Green tensor 
G*,^v(»"A, fA, to) with position va at the center of a small 
spherical cavity of radius i?c which is embedded in an 
infinitely extended bulk material of permittivity ea and 
permeability 11 a- The respective tensor L^?j^(r^, r^i, w) 
reads [H 



Lcil(^A,rA,t^) 



(Bl) 



where 



HAh[^\z) 


zah[^\zo) 


' -h['\zo) 




^J.Ah[^\z) [zoji{zo)]' ~ ji{zo) 


zh'^^\z) 


1 



(B2) 



[z{) = ujRc/c, z — uazq] the primes indicate derivatives 
with respect to zq and z] , with ji (x) and h^i"^ (x) being the 
first-kind spherical Bessel and first-kind spherical Hankel 
functions. 

The local-field correction factors multiplying G in 
Eqs. (|70p and (j7ip are determined by comparing the 
Green tensor Gcav{f,rA,uj) (with at the center of the 
cavity and r at an arbitrary position outside the cavity) 
with the bulk Green tensor Gbuik('', ''a, t^) of an infinite 
homogeneous medium without the cavity. 
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APPENDIX A: INTERMEDIATE STATES AND 
CORRESPONDING DENOMINATORS IN 
EQ. ([57)) 

Here we list the intermediate states contributing to the 
vdW interaction, Eq. (j57p . and the corresponding energy 
denominators (Tab.[l|. 



Gbuik{r,rA,uj) 



{a(— ifcp)l — b{—ikp)epep} (B3) 



with 



a{x) = 1 + X + x"^, b{x) = 3 + 3x + x"^ 



(B4) 



[k = UAUj/c, p = \r ~ VaI, Bp ^ [r - r^)/p]. In pB 
present case, the required tensor Lcav(^, ■'"A, i^) reads |3g | 



,(r,r^,tj) 



D{uj) {a(-ifcp)' - b{-ikp)epep} (B5) 



where 
D{uj) = 



MA " 


h'tHzo) Nji(zo)]' - ji(2o) 


zoh^i\zo) 


n 


PAh^i\z) [zQii^zo)]' ~ ji{zQ) 


'zh[^\z) 





, (B6) 
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Case 






|m> 


Denominator 




(1) 


|fcA,Ofl>|l(i)> 


|0a,0s>|1(2)1{3)) 


!0a,«b)|1(4)) 


Di^ = 
Dib = 


- {ui'X + Ul'){ui' + iu){ljJB + io'), 
-- {UJ% + Uj'){uj' + iu){Lu'B + Uj) 




(2) 


|fcA,OB>|l(i)> 


\kA,lB}m) 


10a,/s)|1(2)) 


D2 = 


{uj\ + Uj'){ij'X + LUb){(^B + 




(3) 


|fcA,Os>|l(l)> 


\kA,lB)m) 


|fcA,0s>|l(2)> 


Di = 






(4) 


|fcA,Os)|l(i)> 


|fcA,«i3)ll{2)l(3)> 


\0a,Ib)\1^4)) 


D4 = 


{u;\ + t^')(^A +UJb+^' + ^)i^^B 


+ u;') 


(5) 


|fcA,Ofl>|l(i)> 


|fcA,«i3)|l{2)l(3)> 


|fcA,0s>|l(4)> 


D5 = 


{Uj\ + Uj'){lo'X +UJb + '^)('^A 




(6) 


\0a,Ib)\1w) 


|0a,0b)|1(2)1(3)> 


|fcA,0fl>|l(4)> 


Dea = 

Deb = 


- {lub + iu'){ijj' + ijj){uja + <^')> 
= {wg + a;')(tj' + Cj)(t^A + 




(7) 


\0a,Ib)\1(i)) 


\kA,lB}m) 


!fcA,0fl>|l(2)> 


Dr = 


{lo'b + (^')i(^A + '^B)i^A + ^) 




(8) 


|OA,iB)|l{l)) 


\kA,lB)m) 


\0a,Ib)\1(2)) 


Ds = 


{lUb + Ijj')(ijJa + iJJB){^B + 




(9) 


|0A,Zi3)|l{l)) 


|fcA,iB)|l{2)l(3)> 


|fcA,0s>|l(4)> 


Dg = 


[Ug + ^'){^A +iOB+i^' + '^)(t^,4 




(10) 


|0A,Zi3)|l{l)) 


|fcA,/s)|l(2)l(3)> 


\0a,Ib)\1(4)) 


Dio = 


= {uj'b + C^')(^A +l^B+^' + 


3+^) 



TABLE I: The intermediate states contributing to the two-atom vdW interaction according to Eq. (I57|l together with the energy 
denominators, where we have used the short-hand notations |1{^)) — \l>.^i^{rtj,,ujii)), |1(^)1(^)) — \l\^i^{rn,uJii) l>.^i^{r,j,ujv)) ■ 



and from Eq. (jB3p . Lbuikl^", ''"A, t^) can be found to be 



P™(x) as 



LbuikCi", rA, uj) = ^^^3 {a{~ikp)\ - b{-ikp)epep} . 

^'^ (B7) 
Comparing Eqs. (|B5|) and (|B7p . we can conclude that, 
on using similar arguments as in Refs. [T^HSl) the mag- 
netic local-field correction factor is given by D / {^An\). 



sm cos 
,(i)n ^dP™(cos6') cos 
da sm 



Combining this with Eg. (IBip and following the line of 
reasoning of Refs. we expand all the terms within 

leading order in uiRc/c to obtain the local-field corrected 
tensors Lioc and in the form of Eqs. (ffi)) and ([75)) . 
Equation ([75]) follows in complete analogy. 



APPENDIX C: GREEN TENSORS L AND K FOR 
A SPHERE 

The free-space part L*-°^ of the magnetic-magnetic 
tensor is the special case e — fi — 1 of the respec- 
tive bulk Green tensor (|B7|) : it obviously coincides with 
— (cj/c)^G [which is a special case of the bulk Green 
tensor (jB3l) ]. According to its definition the scat- 

tering part of L can be found from [3^] 



r-(i)/ ' ^ i^ov^ 2n+l ^ (n-m)! 
>{r,r ,u})^— > ^ — — — > 2-^o„ W -— 

1 ) [n + my. 

X V [Bf (c^)M„™,p(r,A:o)M„„,,p(r',fco) 



47r n(n + 1) -^-^ 

n— 1 m— 

p=±i 

+B^(c^)Af„™^p(r,fco)iVnm,p(r',fco)] , (CI) 

where B^' and are defined by Eqs. and pTU)l . 

are even (p — -1-1) and odd (p — —1) 
spherical wave vector functions and in spherical coor- 
dinates can be expressed in terms of spherical Hankel 
functions of the first kind hn"^ (x) and Legendre functions 



,±i(r,fc) = ^^^^/iW(fcr)i^r(cos0)^°'(m0)e, 



kr 

1 dirh'n^kr)] 



AP^{cose) cos 
d^ sin 



{m4))eg 



7~n m n 

sm cos 
They are related to each other via 

V X M„m^±i(r, k) = kNnrn,±i{r, k), 

V X Ar„„,.,±i(r, fc) = kMn,n.±i{r, k). 



(C3) 



(C4) 
(C5) 



Substituting Eq. (|C1|) into Eq. ([5i)) and making use 
of the relations (jC4p and (jCSp one sees that L^^-* and 
— (a;/c)^G^^^ [and consequently L and —{uj/cYG\ can be 
converted into one another by interchanging i?*^ and 
, or equivalently interchanging e and /i. With this 
knowledge, a comparison between Eqs. ([SO)) and ([55]) re- 
veals that J7mm may be obtained from Uee by replacing 
a with and interchanging e <-> /i. 

The scattering part of tensor K may be found by 
substituting Eq. (jCl[) in ([67)l and making use of rela- 
tions EH) and (IC5l). 



K«(r,r',..) = Mv4!i±L 

^ ' ' ^ 47r ^ n(n + 1) 

n— 1 ^ ^ 



E 

m— 



(2 - (^Om) 



[n + m)\ 



Ar„„,p(r,fco)M„™,p(r',fco) 

p=±i 



(C6) 
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Assuming, without loss of generality, that the coordinate 
system is chosen such that its origin coincides with the 
center of the sphere and the two atoms are located in the 
xz plane as shown in Fig. [TJ 

rA^irA..OA,0), rB^irB..0B,7r), (C7) 



and 



the summations over m and p in Eq. (|C6|) can be per- 
formed in a way similar to Ref . ^3§\ , leading to 



(n— to) 
(n+m) 



-^Nnm..p{rB, ko)Mnm,p{rA,ko) 



m=0 p— ±1 

= -^J^l"-!" + sine P/, (7)6^^60^ 

+ Q,f F„(7)ee,e^^ - Q,f P,;(7)e^,ee^ } , (C8) 



2in(i?w/c)2"+i s{lu) - 1 
[(2n+l)!!]2 £(w) + l' 



(D4) 



and Eqs. (|llip - (|113p approximately reduce to 



In order to illustrate the application of the approxima- 
tion scheme to the tensor K^^-* given by Eq. pOSp let us 
consider, for example, the component -ftT^J,' . Making use 
of Eqs. and jDSj we find that 



E E(2-^o. 



(n — to)! 



Mnm.pirB, ko)Nnm,p{rA,ko) 



m— p— ±1 



(n-|-TO)! 

^ { - QnK h)ees e^^ + n{n + l)Qn sin 6 

(7)60^69^}. (C9) 



Combining Eqs. ([CeD, (|C8j) and (|C9D we arrive at 
Eq. (|108p for the Green tensor. 



X ^i{lu) - 1 



47rP3 ^(^) + 1 (1 - 2tg t2)3/2 

(D7) 



[t — R /(rATB)] where the identity 



oo 



has been used. Recalling condition (|119p . we have 



(D8) 



APPENDIX D: LIMITING CASE OF A LARGE 
SPHERE 

In the limiting case of a large sphere where the condi- 
tions (|119p and (|120p are met, the leading contributions 
to Eqs. (|117p and (jllSp come from terms with n ^ 1 
(see Ref. [32,]), for which the spherical Bessel and Hankel 
functions can be approximated by 



and 



jn{z) = 



(2n-H 1)!! 



.(2n- 1)!! 

„n+l 



1 - 



An + 6 



4n - 2 



(Dl) 



(D2) 



Hence, Eqs. ()109p and (jllOp are approximated by 

M. .^ 2in(P^/c)^"+i ^,{u;)-l 
^« - [(2n+l)!!p ^ ^ 



i*^ = 1 - fc 



-I- 



P2 



i?2 



(D9) 



implying that 

1 - 2i7 + ~ e 



Using Eq. (jPlOp in (jPT)) we end up with 



K'^ArB^rA^uj) - 



= 1^- (DIO) 



47r/:[ ijl{uj) + 1 



(Dll) 



The other components of K^^^ can be evaluated in a sim- 
ilar way. Substituting the resulting expressions for K*-^-* 
into Eqs. (|106p and (jl07p . and summing them in accor- 
dance with Eq. (fTU5)) leads to Eq. ([HT]). 
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